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ABSTRACT 

The  response  of  simply  supported  cross-ply  symmetrically  composite  plates  subjected  to  uniformly  distributed 
load,  with  lamination  [0  90  90  0]  discussed.  Matlab  is  used  to  perform  the  analysis,  depending  on  classical  lamination  plate 
theory.  A  number  of  factors  such  as  aspect  ratio,  side  to  thickness  ratio  and  modulus  ratio  and  their  effect  on  deflection  and 
stresses  of  laminated  composite  plate,  subjected  to  a  uniformly  distributed  load  have  been  studied.  The  results  showed  that, 
the  effect  of  coupling  is  to  increase  the  deflection  to  increase  the  aspect  ratio  and  modulus  ratio,  and  increase  the  stress 
with  an  increase  in  the  side  to  thickness  ratio  and  modulus  ratio. 
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INTRODUCTION 

Composite  materials  are  used  in  different  systems  such  as  space  structures,  automobiles,  sport  equipment  and 
electronic  circuit  boards.  These  materials  are  efficiently  in  applications,  that  required  high  strength  to  weight  and  stiffness 
to  weight  ratios. 

Fiber-  reinforced  composites  are  manufactured  in  laminates.  This  laminate  consists  of  individual  lamina  or  plies. 
There  are  many  factors  that  effect  on  laminate  response  to  applied  loads.  Fiber  material  and  orientation,  fiber  and  matrix 
material  and  layer  sequence  are  some  of  the  variables,  affecting  the  response  of  a  laminate. 

Analysis  of  deformation  in  composite  structure  is  of  fundamental  importance,  in  the  experimental  determination 
of  the  layer  properties  and  exact  solutions  are  useful  in  developing  a  numerical  model. 

Xu  and  Wu  [11  explained  a  two  dimensional  analytical  solution  for  simply  supported  composite  beams  with 
interlayer  slips  by  consideration  of  a  shear  deformation  effect.  Shen  [21  gave  the  nonlinear  analysis  for  bending  of  simply 
supported  functionally  graded  nanocomposite  plates,  subjected  to  transverse  uniform  or  sinusoidal  loads  in  thermal 
environments  and  for  this  purpose,  he  used  a  higher  order  theory  to  derive  governing  equations. 

Torabizadeh  and  Fereidoon  13  presented  an  analytical  and  numerical  solution  for  general  laminated  and  thermal 
loading,  based  on  classical  lamination  plate  theory  (CLPT).  Sidda,  Ramanjancy,  Suresh  and  Vijay  used  a  finite  element 
analysis,  to  study  the  effect  of  transverse  shear  deformation  on  deflection  and  stress  of  laminated  composite  plates. 
Marina  [4]  presented  the  equations  of  bending  in  an  arbitrary  cross  section  of  simply  supported  plates,  under  distributed 
load  based  on  a  partial  layer  wise  theory,  which  is  based  on  the  assumed  displacement  field  in  the  form  of  the  double 
trigonometric  Fourier  s  series. 

Junaid,  Agarwal  and  Vikas  [5]  studied  the  behavior  of  FSDT  plates  under  transverse  loading  condition  and 
estimated  the  influence  of  stacking  sequence,  fiber  orientation,  layer  thickness,  aspect  ratio  and  the  number  of  layers  in  the 
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laminated  composite  plates  on  the  prediction  of  maximum  strength  of  the  plates. 

The  objective  of  this  investigation  is  to  investigate  the  response  of  simply  supported  composite  plates  subjected  to 
uniformly  distributed  load. 

The  Problem  Statement 

A  rectangular  composite  plate  with  length  a,  width  b  and  thickness  h,  consisted  of  4  plies  was  considered. 
Ply  orientation  were  in  0  and  90  degrees  and  the  stacking  sequence  of  the  laminate  was  cross-plied  only  as  in  Figure  (1) 


Figure  1:  Schematic  View  of  Composite  Plate 

THEORETICAL  FORMULATION 

Governing  Equation 

The  governing  equation  for  static  bending  in  the  absence  of  thermal  effect  and  in  plane  forces  is 


[6] 


D^+2tD12+2D66^2  +  D22^=q  1 

Where  D  (Bending  stiffness) 

W  (Displacement  along  the  coordinate-z-) 
q  (The  applied  load) 

The  simply  supported  boundary  conditions  on  all  four  edges  of  the  rectangular  plates  are. 

w  =  0at;r  =  0,a  and  y  =  0  ,  b  2 

Where  a  and  b  are  the  dimensions  of  the  plate  in  x  and  y  respectively. 

Deflection  function 

The  displacement  w  and  the  applied  load  q  (x,  y)  in  the  Navier  method  are  expanded  in  a  double  trigonometric 
(Fourier)  series  of  unknown  parameters.  The  simply  supported  boundary  conditions  for  rectangular  plate  in  equation  (2)  are 
satisfied  by  the  following  forms  of  transverse  deflection  and  applied  load. 


w(x,y)  =  2£=1'Zm=iWmn  sinax  sinfiy 
y)  =  Yln=i  Zm=i  Qmn  sinax  sinfiy 
Where  a  —  mn/a  and  /?  =  nn/b 


3 

4 


Wmn  Are  coefficients  to  be  determined  such  that  the  governing  equation  -1-  is  satisfied  everywhere  in  the  domain 
of  the  plate. 
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Qmn  Are  load  coefficient  and 

Qmn  =  ^  Jo  /;  y)  sinax  sinfiy  dx  dy  5 

Substituting  of  equations  -3-, -4  -and  -5-  in  equation  -1-  results 

£n=i£m=i{-MW0n«4  +  2 (D12  +  2 D66)a4/?4  +  C22/?4]  +  Qmn)  sinax  sinpy=0  6 

The  equation  above  must  be  satisfied  for  every  point  (x,y),  0  <  x  <  a  and  0  <  y  <  b,  and  the  express  inside  the 
brackets  should  be  zero  for  every  m  and  n. 

W  =  - -  7 

mn  D11a'»+2(D12+2D66)a4/?4+022/?4 


Let. 

D11a 4  +  2  (D12  +  2  D66)a4/?4  +  D22/?4  =  dmn  8 

Therefore  Wmn  9 

d-mn 

Then  the  solution  in  equation  -3-  becomes 

w(x,y)  =  EUEU  j^sinaxsinpx  10 

“mn 

For  uniformly  distributed  load  q(x,y)  —  q0  by  instituting  in  equation  -5-  and  integrated  the  Navier  solution  will 
be 


Qmn  =  0  f°r  m,  n,  even 

And  Qmn  —  176i?°  for  m,  n,  odd 
Jr2mn 


Substituting  equation  -12-  in  equation  -9-results 
Wmn  =  716<?° —  for  m,  n  odd 

n  mndmu 


n 

12 


13 


And  therefore  the  displacement  in  simply  supported  composite  plate  under  uniformly  distributed  load  will  be 


w{x,y) 


Y’OD  Y’OO  Ibt Jq 

2m=l  2im=l  Jl2mndmn 


sinax  sin  fix 


For  thin  plates  the  in-plane  stresses  can  be  computed  from  equation  171 
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(Q11Q12  oy 
— z ] Q12Q22  o| 
1  0  0  q-66  ) 


f  dzw 
dx2 
d2w 
dy2 
j  d2w 
V  dxdyJ 
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By  differentiating  w  and  substituting  in  equation -15- yields 
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=  z Vx  V®  W 
Aim=l  Lm= 1  rv  n 


(QT«2  +  «;f^2)Sin=Sin  = 


+  P2)sin 

-(«■) 


tyitzx  .  nny  \ 

- 5m  — 

a  b 


o  n  mnx  nny 

-ZQee  Japcos—  COS  — 


—  (k) 

Where:  Q..  Transformed  reduced  stiffness  matrix. 
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Validation  of  the  Analysis 

In  order  to  verify  or  validate  the  theoretical  analysis  used  here  a  comparison  with  a  numerical  example  of  [81 
considered.  The  boundary  condition  is  simply  supported  and  with  the  following  material  and  geometric  properties 

El/Ez  =  25,  G  =  0.  5  E2,  v12=  0.3,  %  =  1  ,  a/h  =  10 


The  results  are  introduced  here  in  non-dimensional  form  using  the  following. 


E2hs  _  h 

w  =  w  *  — : —  ;  gx  -  ax  — ; 

«4Qo  aqo 


(Tv  —  <TV 


aqQ 


aqQ 


Table  1  represents  a  comparison  of  result  of  non-dimensional  deflection  obtained  from  [81  as  a  finite  element 
analysis  and  from  the  current  analysis. 


Table  1:  Non-Dimensional  Displacement 


Analysis 

w 

FEM  Reference  (8) 

0.00699 

Used  analytical  analysis 

0.0068 

Difference  (%) 

2.7 

RESULTS  AND  DISCUSSIONS 


The  classical  lamination  plate  theory  used  here  for  analysis  of  simply  supported  composite  plate  subjected  to 
uniformly  distributed  load  which  is  with  the  width  (b  =  10  m),  thickness  (h  —  0.1  m)  and  aspect  rai  (a/ E  —  lto  5)  also 

side  to  thickness  ratio  (a/fl  =  10  to  40  ).  The  applied  load  q0  was  (1  ^/m2  )•  'n  this  analysis  4  plies  (0  90  90  0)° 

symmetric  composite  (Reddy)  [9]  considered  with  fowling  material  properties  (E2  =  1  *  106  ,  v12  =  0.25,  and 

C12  =  0.  5  *  E2)  and  the  other  properties  assumed  to  be  as  (  E3=E2  C13  =  C12,  v23  =  v13=  v12  *  —  )  and  different 

'  E  i 

Eii 

modulus  ratios  (  1/ =  1  to  35  ).  The  results  obtained  for  non  dimensional  deflections  (w  ),  normal  stresses  (ox  ,  ay ) 
and  shear  stresses  ( r~y )  are  plotted  as  a  function  of  aspect  ratio  (a/^)  For  different  modulus  ratios  as  in  Figures  2,  3,  4  and 
5  respectively.  These  figures  show  the  effect  of  bending-stretching  coupling  and  plate  aspect  ratio  on  deflections  and 
stresses. 

From  Figure  2,  it  is  obvious  that  non  dimensional  deflection  (vv_)  is  the  maximum  at  an  aspect  ratio  (a/^  =  1), 

E  E 

modulus  ratio  (  1/ ^  =  1)  and  minimum  at  aspect  ratio  (a/j,  =5),  modulus  ratio  (  1/E  —  35  ).  This  behavior  is  related 

£ 

to  young’s  modulus,  where  increasing  the  modulus  ratio  Y^causes  the  coupling  coefficient  to  increase  and  this  causes 
the  nondimensionl  deflection  (w-)  to  decrease.  The  effect  of  coupling  is  signed  for  a/^less  than  about  3.5  and 
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insignificant  for  all  values  greater  than  3.5. 

From  Figures  3-5  can  observe  the  effect  of  coupling  is  to  decrease  the  stress  with  increasing  the  aspect  ratio, 

The  normal  nondimensional  stresses  a~  ,  <7y  and  shear  non  dimensional  stress  r“yare  maximum  at 

E  E 

,  1/g  —  35and  minimum  ata/^  =  5,  X  / E2  ~  ^  Fh i&  is  because  the  area  of  the  composite  plate  increases  as  the  aspect 
ratio  increase,  therefore  the  applied  load  per  unit  area  is  decreased. 

Figures  6-9  show  the  effect  of  modulus  ratios  and  the  side  to  thickness  ratios  on  deflection,  normal  stresses  and 
shear  stresses  for  simply  supported  composite  plates  laminate  [0  90  90  0],  under  uniform  distributed  load. 
The  side  to  thickness  ratio  has  less  effect  on  the  deflection,  especially  for  larger  modulus  ratios,  but  has  a  considerable 
effect  on  stresses.  The  deflection  decreases  with  increasing  the  modulus  ratio  and  side  to  thickness  ratio,  while  the  stresses 
decrease  with  the  modulus  ratios  and  increase  with  increasing  the  side  in  thickness. 

£ 

The  nondimensional  deflection  is  maximum  for  side  to  thickness  ratio  (a/^  =  35),  modulus  ratio  (  X  / e2  =  5), 

and  minimum  of  side  to  thickness  ratio  (a/^  =  10),  modulus  ratio  (  1/g  =  53)  and  this  can  reduce  with  increasing  the 

plate  area  with  increasing  side  to  thickness  ratio,  therefore,  that  can  causes  the  applied  load  per  unit  area  decreases. 
The  normal  stress  (rr“)  is  maximum  at  side  to  thickness  ratio  (a/^  =  35)  without  considerable  effect  of  the  modulus  ratio, 

while  ( <Ty )  and  (rjy)  is  maximum  at  side  to  thickness  ratio  (a/^  =  35),  modulus  ratio  (  1 /g  —  5)  and  minimum  of  side 

£ 

to  thickness  ratio  (a/fl  —  10),  modulus  ratio  (  1  / g  —  35). 


Aspect  ratio, (a/b) 


Figure  2:  Nondimensionalized  Maximum  Transverse  Deflection 
(W)  Versus  Aspect  Ratio  (a/b)  For  Different  Modulus  Ratios  E1/E2 
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Figure  3:  Nondimensionalized  Normal  Stress  (ex')  Versus  Aspect  Ratio 
(a/b)  for  Different  Modulus  Ratios  E1/E2 
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Aspect  ratio, (a/b) 


Figure  4:  Nondimensionalized  Normal  Stress  (ay')  Versus  Aspect 
Ratio  (a/b)  for  Different  Modulus  Ratios  E1/E2 


Figure  5:  Nondimensionalized  Shear  Stress  (r  '„)  Versus  Aspect  Ratio 
(a/b)  For  Different  Modulus  Ratios  E1/E2 


Figure  6:  Nondimensionalized  Deflection  (W")  Versus  Modulus  Ratio 
(Ei/E2)  For  Different  Side  to  Thickness  Ratios  a/h 
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Figure  7:  Nondimensionalized  Normal  Stress  (a  ',)  Versus  Modulus  Ratio 
(Ej/Ej)  For  Different  Side  to  Thickness  Ratios  a/h 


Figure  8:  Nondimensionalized  Normal  Stress  (a  ',)  Versus  Modulus  Ratio 
(E^Ej)  For  Different  Side  to  Thickness  Ratios  a/h 


Modulus  ratio, (E^Eg) 

Figure  9:  Nondimensional  Shear  Stress  (rxy)  Versus  Modulus  Ratio 
(Ei/E2)  For  Different  Side  to  Thickness  Ratios  a/h 

CONCLUSIONS 


By  using  a  classical  lamination  plate  theory  and  mat  lab  code  at  various  aspect  ratios,  side  to  thickness  ratios  and 
modulus  ratios  were  performed  to  investigate  how  it  affects  the  deflections  and  stresses.  The  results  are  plotted  in  previous 
figures  and  observed  that,  the  deflections  are  larger  for  small  aspect  ratios  and  modulus  ratios;  the  side  to  thickness  ratios 

£ 

has  an  effect  on  the  deflection  for  larger  ratios  of  1/g  . 

The  deflection  decrease  and  the  stress  increase  with  increase  of  modulus  ratios  and  side  to  thickness  ratios. 
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